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ABSTRACT 

Context. The frequency of brown-dwarf companions in close orbit around Sun-like stars is low compared to the frequency of plane- 
tary and stellar companions. Currently, there is no comprehensive explanation for this lack of brown-dwarf companions. 
Aims. Through the combination of the orbital solution obtained from a stellar radial-velocity curve and Hipparcos astrometric mea- 
surements, we attempt to determine the orbit inclination and therefore the mass of the orbiting companion. By determining the masses 
of potential brown-dwarf companions, we can refine the properties of the companion mass-function. 

Methods. The radial-velocity solutions revealing potential brown-dwarf companions are obtained for stars from the CORALIE and 
HARPS planet-search surveys or from the literature. The best Keplerian fit to our radial-velocity measurements is found by adjust- 
ment with the Levenberg-Marquardt algorithm. The spectroscopic elements of the radial-velocity solution constrain the fit to the 
Intermediate Astrometric Data of the new Hipparcos reduction. The astrometric solution and the orbit inclination are found using 
non-linear ^-minimisation on a two-parameter search grid. The statistical confidence of the adopted orbital solution is evaluated 
based on the distribution-free permutation test. 

Results. The discovery of 9 new brown-dwarf candidates orbiting stars in the CORALIE and HARPS radial-velocity surveys is re- 
ported. New CORALIE radial velocities yielding accurate orbits of 6 previously-known hosts of potential brown-dwarf companions 
are presented. Including the literature targets, 33 hosts of potential brown-dwarf companions are examined. Employing innovative 
methods, we use the new reduction of the Hipparcos data to fully characterise the astrometric orbits of 6 objects, revealing M-dwarf 
companions with masses between 90 Mj and 0.52 M B . Additionally, the masses of two companions can be restricted to the stellar 
domain. The companion to HD 137510 is found to be a brown dwarf. At 95 % confidence, the companion of HD 190228 is also a 
brown dwarf. The remaining 23 companions persist as brown-dwarf candidates. Based on the CORALIE planet-search sample, we 
obtain an upper limit of 0.6 % for the frequency of brown-dwarf companions around Sun-like stars. We find that the companion-mass 
distribution function is rising at the lower end of the brown-dwarf mass range, suggesting that in fact we are detecting the high-mass 
tail of the planetary distribution. 

Conclusions. Our findings agree with the results of previous similar studies and confirm the pronounced paucity of brown-dwarf 
companions around Sun-like stars. They are affected by the Hipparcos astrometric precision and mission duration, which limits the 
minimum detectable companion mass, and some of the remaining candidates probably are brown-dwarf companions. 
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1. Introduction 

The mass distribution of close companions to Sun-like stars ex- 
hibits a clear deficit of brown dwarfs compared to planets and 
stellar binaries (Grether & Lineweaver 2006). With a frequency 
of less than 1 % (Marcy & Butler 2000), brown-dwarf compan- 
ions to solar-type stars at separations below 10 AU are much less 
common than planetary companions with a frequency of ~1 % 
(e.g. Udry & Santos 2007) and stellar binaries with a frequency 
of -13 % (Duquennoy & Mayor 1991; Halbwachs et al. 2003) 

* Based on observations collected with the CORALIE Echelle spec- 
trograph on the Swiss telescope at the European Southern Observatory 
in La Silla, Chile, on observations made with the ESA Hipparcos as- 
trometry satellite, and on observations made with the HARPS instru- 
ment on the ESO 3.6-m telescope (GTO programme 072.C-0488). 
The CORALIE and HARPS radial-velocity measurements discussed in 
this paper are only available in electronic format at CDS via http: 
//cdsweb.u-strasbg . fr/cgi-bin/qcat?J/A+A/. 



at similar separations. Brown dwarfs are substellar objects, mas- 
sive enough to burn deuterium but too light to permit hydro- 
gen burning, which places them in the mass range of approxi- 
mately 13-80 Jupiter masses (Burrows et al. 1997; Chabrier 
& Baraffe 2000). They therefore constitute objects at the tran- 
sition between planets and stars, though the lower end of the 
brown-dwarf mass range overlaps with the one of massive plan- 
ets (e.g. Luhman et al. 2009), and the distinction between planets 
and brown dwarfs may require to trace the individual formation 
process. A comprehensive explanation for the formation of low- 
mass objects, covering the range of planets to very low-mass 
stars, remains to be found, although the progress in this direc- 
tion is formidable as reported in the reviews of Whitworth et al. 
(2007), Luhman et al. (2007) and Burgasser et al. (2007). 

Recently, a few companions with masses of 18 - 60 Jupiter 
masses (Mj) were determined through astrometry (e.g. Martioli 
et al. 2010; Benedict et al. 2010) or discovered in transiting 
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systems (e.g. Deleuil et al. 2008, Anderson et al., in prepara- 
tion, Bouchy et al., in preparation). Yet, most potential brown 
dwarf companions are discovered in radial velocity surveys (e.g. 
Nidever et al. 2002). In these high-precision programs, aiming 
at planet detection, they are easily spotted and characterised be- 
cause of their large signatures, by far exceeding the radial veloc- 
ity measurement precision (e.g. Patel et al. 2007; Bouchy et al. 
2009). However, radial-velocity measurements alone do not con- 
strain the orbit inclination. Therefore, they can not reveal the 
companion mass, but yield a lower limit to it. Observations of 
complementary effects like the transit light-curve (Deleuil et al. 
2008) or the astrometric motion of the host star are required to 
solve this ambiguity and to determine the companion mass. 

High-precision astrometry with the Hubble space telescope 
(HST) has been used to find the masses of planetary compan- 
ions discovered by radial velocity (Benedict et al. 2010, and 
references therein) and has demonstrated the power of astrom- 
etry to distinguish between planetary, brown-dwarf, or stel- 
lar companions. Ground-based optical interferometers achieve 
the necessary precision to detect orbital motions of binary and 
multiple systems and to solve for the orbital parameters to- 
gether with radial-velocity measurements (Hummel et al. 1993; 
Muterspaugh et al. 2006; Lane et al. 2007). The astromet- 
ric detection of planetary orbits using IR-interferometry is the 
aim of the extrasolar-planet search with PRIMA project (ESPRI, 
Launhardt et al. 2008). On a larger scale, the GAIA satellite (e.g. 
Lindegren 2010) will advance the field of astrometric detection 
and characterisation of planetary systems, due to its outstanding 
measurement precision (Casertano et al. 2008). 

Hipparcos astrometry has been extensively used to constrain 
the masses of sub-stellar and stellar companions in multiple sys- 
tems (Perryman et al. 1996; Pourbaix & Jorissen 2000; Zucker 
& Mazeh 2001; Torres 2007). Han et al. (2001) used Hipparcos 
astrometry to analyse 30 extrasolar-planet candidates and con- 
cluded that half of these companions were rather brown of M 
dwarfs. With the consequent rectification by Pourbaix & Arenou 
(2001), showing that the results of Han et al. (2001) had been 
biased by the fitting procedure, it became clear that the fitting of 
astrometric data has to be done very carefully, especially when 
the size of the orbital signature approaches the measurement pre- 
cision of the instrument (Pourbaix 2001). All of these past works 
used the original Hipparcos data (Perryman et al. 1997). 

In 2007, the new reduction of the raw Hipparcos data (F. van 
Leeuwen 2007) has been released, which represents a signifi- 
cant improvement in quality compared to the original data (van 
Leeuwen 2007). It also includes a revised version of the interme- 
diate astrometric data in a new format that facilitates the search 
for the signature of orbital motion. We use the new Hipparcos re- 
duction for our study, which aims at determining the masses of 
companions to stars, that were detected by high-precision radial 
velocity measurements and found to have minimum masses in 
the brown-dwarf domain. Our target sample is composed of stars 
from the CORALIE and HARPS planet surveys. In addition, we 
included a list of targets selected from the literature. To validate 
our method of astrometric orbit determination, we also analysed 
a few comparison targets with published astrometric orbits. This 
study further explores the brown dwarf desert and is in line with 
the foregoing works by Halbwachs et al. (2000) and Zucker & 
Mazeh (2001). 

The paper is organised as follows. The brown-dwarf candi- 
dates from the CORALIE and HARPS surveys are presented in 
Sect. 2 together with 18 candidates selected from the literature. 
The method of combining radial-velocity orbits and Hipparcos 



Table 1. Hipparcos parameters of the surveyed stars 



Nr. 


Object 


Sp. T. 


V 


B-V 




tzr 




HD 


HIP 








(mas) 


1 


3277 


2790 


G8V 


7.59 


0.73 


34.7 


±0.7 


2 


4747 


3850 


G8 


7.30 


0.77 


53.5 


±0.5 


3 


17289 


12726 


GOV 


7.56 


0.59 


20.8 


± 1.2 


4 


30501 


22122 


K0V 


7.73 


0.88 


47.9 


±0.5 


5 


43848 


29804 


K1IV 


8.80 


0.93 


26.4 


±0.8 


6 


52756 


34052 


K1V 


8.80 


1.18 


57.4 


± 1.2 


7 


53680 


33736 


K5V 


8.61 


0.90 


30.1 


±0.6 


8 


74014 


42634 


K0 


7.73 


0.76 


29.6 


±0.7 


9 


89707 


50671 


G1V 


7.29 


0.55 


30.6 


±0.6 


10 


154697 


83770 


G5 


7.97 


0.73 


30.6 


± 1.0 


11 


164427 A 


88531 


G4IV 


7.01 


0.62 


26.3 


±0.7 


12 


167665 


89620 


F8V 


6.48 


0.54 


32.8 


±0.5 


13 


189310 


99634 


K2V 


8.61 


0.90 


31.0 


±0.6 


14 


211847 


110340 


G5V 


8.78 


0.66 


19.8 


± 1.3 


15 




103019 


K5 


10.45 


1.33 


28.3 


±2.3 



astrometry is described in Sect. 3. The results are summarised in 
Sect. 4 and discussed in Sect. 5. We conclude in Sect. 6. 

2. Brown dwarf candidates 

Our target sample consists of 33 stars, that exhibit radial velocity 
variations caused by a companion with minimum mass Mi sin i 
in the brown-dwarf mass range of 13 - 80 My. It contains 14 
stars from the CORALIE survey (Udry et al. 2000) and one star 
from the HARPS planet search program (Mayor et al. 2003). To 
date, the CORALIE planet survey has contributed to the discov- 
ery of more than 50 extrasolar planets (e.g. Naef et al. 2001; 
Mayor et al. 2004; Tamuz et al. 2008). CORALIE is an optical 
echelle spectrograph mounted on the 1.2 m Swiss Telescope lo- 
cated at the European Southern Observatory in La Silla, Chile. 
A description of the instrument can be found in Queloz et al. 
(2000) and the references therein. With its unmatched preci- 
sion, the HARPS instrument permitted the discovery of more than 
70 planetary companions to date, including earth-mass planets 
(e.g. Pepe et al. 2007; Santos et al. 2007; Mayor et al. 2009; 
Forveille et al. 2009). The characteristics, radial velocities, and 
orbital solutions of the 15 CORALIE and HARPS stars are pre- 
sented and discussed in Sects. 2.1-2.3. Additionally, 18 poten- 
tial brown-dwarf host stars with radial -velocity orbits were se- 
lected from the literature and are listed in Sect. 2.4. These in- 
clude HD 190228, initially a planet-host star (Sivan et al. 2004), 
whose importance we recognised during the analysis. 

2.1. Stellar characteristics 

The identifiers and basic stellar characteristics of the CORALIE 
stars and of HIP 103019 from the HARPS survey are listed in 
Table 1. The apparent visual magnitude V, the colour B-V, and 
the parallax vj are from the new Hipparcos reduction (F. van 
Leeuwen 2007), whereas the spectral type is from the original 
Hipparcos catalogue (Perryman et al. 1997). Table 2 displays the 
derived stellar parameters. The uncertainty of V is below 0.002 
mag. The absolute magnitude My is derived from the Hipparcos 
magnitude and parallax. The effective temperature T e g , the sur- 
face gravity logg, and the metallicity [Fe/H] are derived from 
the spectroscopic analysis of high-signal-to-noise spectra with 
the method presented in Santos et al. (2004) and using the Fel 
and Fe II lines listed in Sousa et al. (2008). The stellar rotation 
parameter v sin i is derived from the calibration of the CORALIE 
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Table 2. Derived stellar parameters of the surveyed stars 



Object 


Instr* 


M v 




Teff 

(K) 




log.? 
(cgs) 




[Fe/H] 

(dex) 


v sin i 
(kms" 1 ) 




Mi 

(M a ) 




Age 
(Gyr) 


HD 3277 


C 


5.30 


5539 


±50 


4.36 


±0.10 


-0.06 


±0.08 


< 2.0 


0.91 


±0.03 


0.4 


- 8.0 


HD 4747 


C 


5.94 


5316 ± 


50" 


4.48 ± 


0.10" 


-0.21 ± 


: 0.05" 


0.79 ± 0.06" 


0.81 


±0.02 


0.1 


-7.3 


HD 17289 


C 


4.15 


5924 


±50 


4.37 


±0.10 


-0.11 


±0.06 


< 2.0 


1.01 


±0.03 


6.0 


-9.3 


HD 30501 


C 


6.13 


5223 


±50 


4.56 


±0.10 


-0.06 


±0.06 


< 2.0 


0.81 


±0.02 


0.8 


-7.0 


HD 43848 


C 


5.91 


5334 


±92 


4.56 


±0.15 


0.22 


±0.06 


< 2.0 


0.89 


±0.02 


0.0 


-5.1 


HD 52756 


C 


7.60 


5216 


±65 


4.47 


±0.11 


0.13 


±0.05 


< 2.0 


0.83 


±0.01 


0.2 


- 2.7 


HD 53680 


C 


6.00 


5167 


±94 


5.37 


±0.29 


-0.29 


±0.08 


2.08 ±0.31 


0.79 


±0.02 


0.7 


-9.4 


HD 74014 


C 


5.08 


5662 


±55 


4.39 


±0.10 


0.26 


±0.08 


< 2.0 


1.00 


±0.03 


0.1 


-3.8 


HD 89707 


C 


4.71 


6047 


±50 


4.52 


±0.10 


-0.33 


±0.06 


< 2.0 


0.96 


±0.04 


0.5 


-6.0 


HD 154697 


C 


5.40 


5648 


±50 


4.42 


±0.10 


0.13 


±0.06 


< 2.0 


0.96 


±0.02 


0.0 


-3.3 


HD 164427 A 


C 


4.11 


6003 


±50 


4.35 


±0.10 


0.19 


±0.06 


< 2.0 


1.15 


±0.03 


1.2 


-4.1 


HD 167665 


C 


4.06 


6224 


±50 


4.44 


±0.10 


-0.05 


±0.06 


< 2.0 


1.14 


±0.03 


0.7 


-3.6 


HD 189310 


C 


6.07 


5188 


±50 


4.49 


±0.10 


-0.01 


±0.06 


< 2.0 


0.83 


±0.02 


0.0 


-7.5 


HD 211847 


C 


5.26 


5715 


±50 


4.49 


±0.10 


-0.08 


±0.06 


< 2.0 


0.94 


±0.04 


0.1 


-6.0 


HIP 103019 


H 


7.61 


4913 + 


: 115 


4.45 


±0.28 


-0.30 


±0.06 




0.70 


±0.01 


0.0 


-7.9 



Notes. (fl) Data from Santos et al. (2005). w C and H stand for CORALIE and HARPS, respectively. 



or HARPS cross correlation function by Santos et al. (2002). 
Finally, the stellar mass of the primary Mi and the age are es- 
timated from the theoretical isochrones of Girardi et al. (2000) 
and a Bayesian estimation method described in da Silva et al. 
(2006) 1 . Stellar ages of main-sequence dwarfs are usually not 
well constrained and we quote the 1-cr confidence interval. The 
spectroscopic analysis was problematic for HIP 103019, which 
may be due to its faintness and late spectral type. Errors on the 
parameters of this star are possibly underestimated. The stellar 
characteristics of the targets from the literature can be found in 
the respective references given in Sect. 2.4. 

2.2. Radial-velocity measurements and orbital solutions 

Optical high-resolution spectra of the stars in Table 2 were 
collected with the CORALIE and HARPS spectrographs over 
time spans extending to 11 years. Radial velocities are esti- 
mated from the cross-correlations of the extracted stellar spectra 
with numerical templates, which depend on the target's spec- 
tral type. During CORALIE observations, a reference Thorium- 
Argon spectrum is recorded simultaneously with the stellar spec- 
trum and is used to measure and correct for residual zero-point 
drifts (Baranne et al. 1996). Photon noise limits the obtained 
precision to typically 2-4 ms _I per epoch. To account for sys- 
tematic drifts, an external error of 5 ms~' is quadratically added 
to CORALIE radial-velocity uncertainties before performing the 
period search and the model adjustment. The abbreviations C98 
and C07 refer to the CORALIE instrument before and after its 
upgrade in 2007, respectively (Segransan et al. 2010). The two 
parameters yc98 and ycoi are introduced to account for the possi- 
bly differing velocity offsets of these instrument configurations. 

Because of the stability of the HARPS instrument and the tar- 
geted precision of 2 ms _1 per epoch for the volume -limited pro- 
gramme, HIP 103019 is observed without Thorium-Argon ref- 
erence (cf. Naef et al. 2010). To account for possible drifts, an 
additional error of 0.5 ms 1 is quadratically added to the HARPS 
radial-velocity uncertainties. 

The Keplerian-orbit solution is found by imposing the model 
function Eq. 1 1 with 6 parameters (Ki , e, to, T, P, and y) to the 
radial velocities of a given star. The fit determines the values of 

1 The web interface is http://stev.oapd.inaf.it/cgi-bin/ 
par am 



the period P, the eccentricity <?, the time of periastron passage To, 
the longitude of periastron co, the radial-velocity semi-amplitude 
K\, and the systemic velocity y. In general, measurements with 
uncertainties larger than 15 ms are not used for the model ad- 
justments. 

A single-companion model fits well the data of 1 1 stars with 
reduced-chi values of Xr = 0.7 - 1.3. For 4 stars, HD 17289, 
HD 89707, HD 211847, and HIP 103019, the fit is not as good 
with^ r = 1.7-2.0 and produces residuals of unusually large 
amplitude, without contesting the presence of a large orbital mo- 
tion. Introducing additional model parameters, corresponding to 
a linear or quadratic drift or a second companion does not sig- 
nificantly improve the fit. Therefore, the moderate fit quality has 
to be attributed to either underestimated radial velocity errors or 
to additional noise sources, such as stellar activity. This is dis- 
cussed for the individual targets in the next Section. 

The orbital periods of HD 4747 and HD 211847 are very 
long and not covered by the measurement time-span. Hence, the 
best-fit period is poorly constrained and realistic confidence in- 
tervals have to be derived through a different method. We chose 
a criterion based on the fit-quality metric Xr and its error cr x 
obtained for the best fit. During a second fit, we fix the period 
during the adjustment, such that the resulting ^-value increases 
to Xr + In this way, we obtain upper and lower limits for 
the orbital parameters. They are listed at the bottom of Tables 
3 and 4. The method to search for the astrometric signature 
in the Hipparcos data presented in Sect 3 relies on the precise 
knowledge of the radial-velocity parameters. Because the peri- 
ods of HD 4747 and HD 211 847 are poorly constrained and ad- 
ditionally the orbital coverage of Hipparcos is very low, we do 
not perform the combined astrometric analysis for these targets. 
They are presented as part of the brown-dwarf candidates char- 
acterised with CORALIE. 

Table 3 displays the orbital elements of the best fit solutions 
and Table 4 lists the basic statistics and the derived orbit param- 
eters (see also Fig. 1). Tq is given as reduced Barycentric Julian 
Date (BJD), that is JD* = BJD - 2400000. N RV is the total 
number of measurements, Ar is the observation time-span, S/N 
is the median signal-to-noise ratio per pixel at 550 nm, <cr RV > is 
the mean measurement uncertainty, Xr is the reduced-chi value 
of the fit, G.o.F. is the goodness of fit, <x (0 _ C ) gives the root- 
mean-square dispersion of the fit residuals, ci\ sin / is the min- 
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Fig. 1. Orbital characteristics of the 15 brown-dwarf candidates 
observed with CORALIE and HARPS. Filled symbols indicate the 
companions with stellar mass identified through this study. The 
error bars are often smaller than the symbol size, though the or- 
bits of HD 4747 and HD 21 1847 are not covered and therefore 
less constrained, see the text for details. Top: Minimum compan- 
ion mass and separation. Bottom: Eccentricity and orbital period. 



imum astrometric semimajor axis expressed in unit of length, 
f(m) is the mass function, a le \ is the semimajor axis of the rela- 
tive orbit, and finally Mi sin i is the minimum companion mass. 
The quoted errors correspond to l-<x confidence intervals derived 
from 5000 Monte-Carlo simulations. Figures 2, 3, and 4 show 
the phase-folded radial velocities and the adopted orbital solu- 
tion. All measurements are available in electronic form at CDS 
as indicated on the title page. 



2.3. Details on the surveyed stars 

- HD 3277 (HIP 2790) is discovered to host a potential brown- 
dwarf companion. The star was listed as an uncertain spec- 
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Fig. 4. Radial velocities (top) and residuals (bottom) to the best 
fit, represented by the solid line, of HD 4747 measured with 
HIRES (orange symbols), CORALIE C98 (red symbols), and 
CORALIE C07 (blue symbols). 
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CORALIE 




2600 2800 3000 3200 3400 3600 3800 4000 4200 4400 4600 4800 5000 5200 
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Fig. 5. Radial velocities (top) and residuals (bottom) to the best 
fit, represented by the solid line, of HD 211847 measured with 
CORALIE C98 (red circles), and CORALIE C07 (blue squares). 



troscopic binary in Tokovinin et al. (2006), who did not have 
sufficient radial-velocity measurements to constrain the or- 
bit. 

HD 4747 (HIP 3850) has a potential brown-dwarf compan- 
ion discovered with radial velocities by Nidever et al. (2002). 
We combine CORALIE observations with the KECK HIRES 
radial velocities of Nidever et al. (2002) retrieved via the 
SB9-catalogue (Pourbaix et al. 2004). However, the orbital 
period is very long and not covered by the measurements, 
see Fig. 4. Tables 3 and 4 show the best-fit solution and 
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Fig. 2. Phase-folded radial velocities of 8 stars with potential brown-dwarf companions. Red circles and blue squares indicate 
measurements with CORALIE C98 and C07, respectively. For HD 43848, the orange diamonds show the MIKE measurements. The 
solid lines correspond to the best-fit solutions. The error bars are smaller than the symbol size. 
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Fig. 3. Phase-folded radial velocities of 5 stars with potential brown-dwarf companions. Error bars are smaller than the symbol size. 



the solutions with fixed period such that the ^-value in- 
creases to Xr + 3o~x ~ 1.56. The resulting period range 
is P = 7900 - 29000 days with corresponding minimum 
companion-mass Mi sin / = 39.6 - 58.1 My. 
- HD 17289 (HIP 12726) is discovered to host a potential 
brown-dwarf companion. Goldin & Makarov (2007) inde- 
pendently found an orbital signature from Hipparcos astrom- 
etry alone with a period of 536 + 12 days, which is close to 
our solution (562.1 + 0.3 days). However, the orbital param- 
eters derived by Goldin & Makarov (2007) are less precise 
and less accurate than the radial-velocity solution. 
The dispersion of the fit residuals is abnormally large for 
CORALIE. The detailed analysis of the cross-correlation 



function (CCF) yields that its bisector velocity span and full- 
width-at-half-maximum (FWHM) are correlated with the or- 
bital phase (Fig. 6). As we will show in Sect. 4, the compan- 
ion of HD 17289 is a star with mass of 0.52 M Q . The two 
stars are very close (< 50 mas, Sect. 4.1) and their estimated 
intensity ratio in the visible is about 1:180. Therefore, the 
companion's light is equally picked up by the science fibre 
of CORALIE and its signature appears in the collected spec- 
tra. The spectral lines of the companion distort the CCF de- 
pending on the orbital phase. As expected the CCF width is 
minimal when the radial-velocity curves of both components 
cross and equal the systemic velocity. The detailed modelling 
of this double-lined spectroscopic binary is outside the scope 
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Table 3. Single-companion Keplerian orbital solutions of the surveyed stars 
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of this paper, but we have identified the binary as the cause 
of a phase-dependent radial-velocity bias, which explains the 
excess noise in our measurements. 

- HD 30501 (HIP 22122) is discovered to host a potential 
brown-dwarf companion. 

- HD 43848 (HIP 29804) has a potential brown-dwarf com- 
panion discovered with radial velocity by Minniti et al. 
(2009). An astrometric analysis using Hipparcos and the so- 
lution of Minniti et al. (2009) was performed by Sozzetti 
& Desidera (2010). We combined the published measure- 
ments 2 with the CORALIE velocities to adjust the Keplerian 
model. While the orbital period is well defined, the distri- 
butions of e, K\, To, and u> obtained from the Monte-Carlo 
simulations are bimodal. This is because the lower velocity 
turnover is not well sampled by the measurements. The fit- 
quality in terms of^- 2 is comparable for all Monte-Carlo sim- 
ulations and can not be used to identify the correct solution. 
In Fig. 7 we show the eccentricity and semi-amplitude dis- 
tributions. Because the realisation of the lower-eccentricity 
orbit is more probable (77 %) than the high-eccentricity or- 
bit (23 %), we consider only the orbits with e < 0.78 to 
find the final solution and to perform the astrometric anal- 
ysis. For completeness, we also list the less probable and 
higher-eccentricity solution in Tables 3 and 4. 

- HD 52756 (HIP 33736) is discovered to host a potential 
brown-dwarf companion. 

2 By comparing the published radial velocities to our CORALIE 
measurements, we discovered that the published lulian dates corre- 
sponding to the radial velocities of Minniti et al. (2009) are misprinted 
and offset by -2000 days. 



- HD 53680 (HIP 34052) is discovered to host a potential 
brown-dwarf companion. 

- HD 74014 (HIP 42634) has a potential brown-dwarf com- 
panion discovered with radial velocities by Patel et al. 
(2007). The original orbit is compatible with our solution ob- 
tained from 118 CORALIE velocities covering the complete 
orbit. 

- HD 89707 (HIP 50671) has a potential brown-dwarf com- 
panion announced by Duquennoy & Mayor (1991) and anal- 
ysed by Halbwachs et al. (2000) based on radial veloci- 
ties collected with the CORAVEL and ELODIE instruments. 
The original orbit is compatible with our solution obtained 
from 37 CORALIE velocities. A combined astrometric anal- 
ysis using Hipparcos was performed by Zucker & Mazeh 
(2001). We find slight evidence for re-emission at the bot- 
tom of the Call H absorption line (Fig. 8) and derive an ac- 
tivity index of log R' HK = -4.66 ± 0.03 from a high signal- 
to-noise CORALIE spectrum. Using the relation of Santos 
et al. (2000), we expect an additional activity-induced radial- 
velocity noise of 12 ms -1 , which explains the unusually- 
large fit residuals encountered for HD 89707. Accounting 
for the additional noise during the adjustment yields an ac- 
ceptable fit quality with^ r = 0.91 ± 0.13. 

- HD 154697 (HIP 83770) is discovered to host a potential 
brown-dwarf companion. 

- HD 164427A (HIP 88531)has apotential brown-dwarf com- 
panion announced by Tinney et al. (2001), whose original or- 
bit is compatible with our solution. A combined astrometric 
analysis using Hipparcos was performed by Zucker & Mazeh 
(2001). 
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Table 4. Companion parameters for the surveyed stars 
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Fig. 7. Eccentricities (left) and velocity semi- 
amplitudes (right) of 10000 Monte-Carlo sim- 
ulations performed to find the radial-velocity 
orbit of HD 43848. The bi-modal structure is 
clearly seen in both histograms. The final so- 
lution is found by considering only the lower 
eccentricity realisations, which is done by sep- 
arating the distributions at K[ = 910 ms~' and 
correspondingly at e ~ 0.78. The probability 
to obtain a high- or low-eccentricity solution is 
calculated by counting the number of realisa- 
tions above and below this threshold value and 
amounts to 77 % and 23 %, respectively. 



- HD 189310 (HIP 99634) is discovered to host a potential 
brown-dwarf companion in a short -period orbit of 14.2 days. 

- HD 167665 (HIP 89620) has a potential brown-dwarf com- 
panion announced by Patel et al. (2007), whose original orbit 
is compatible with our solution. 

- HD 211847 (HIP 110340) is discovered to host a poten- 
tial brown dwarf companion. The long-period orbit is not 
completely covered by our measurements (Fig. 5) and the 
CORALIE upgrade occured at minimum radial velocity. To 
obtain a reasonable solution, we therefore fix the offset be- 
tween the parameters yc98 an d Jcoi t0 zero, which is the 
mean offset value found for three other stars of same spectral 
type. As for HD 4747, we derive realistic confidence inter- 
vals by fixing the period during the adjustment, such that 



the resulting ^-value increases to Xi + ~ 2.32. The 
period lower limit is P = 3750 days and we set the up- 
per limit to P = 100000 days, because we cannot derive 
an upper limit using the defined criterion. The correspond- 
ing minimum companion mass is M2 = 17.0 - 23.3 Mj. 
The dispersion of the best-fit residuals is abnormally large 
for CORALIE. We find evidence for re-emission at the bot- 
tom of Ca II H and Ca II K, see Fig. 9, and derive an activ- 
ity index of logR' HK = -4.74 + 0.04 from a high signal- 
to-noise CORALIE spectrum. Using the relation of Santos 
et al. (2000), we expect an additional activity-induced radial- 
velocity noise of 11ms -1 , which explains the fit residuals 
encountered for HD 211847. Accounting for the additional 
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Fig. 6. Bisector velocity span (top) and FWHM of the cross- 
correlation-function {bottom) as function of orbital phase for 
HD 17289. The corresponding radial velocities are shown in 
Fig. 2. Clearly, both observables show a phase-dependent mod- 
ulation, which is likely caused by the companion. 
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Fig. 8. Ca II H emission region for HD 89707 showing slight re- 
emission at the bottom of this absorption line. 



noise during the adjustment yields an acceptable fit quality 
wither = 1.15 ±0.15. 
- HIP 103019 is discovered to host a potential brown-dwarf 
companion using HARPS. The star is included in the search 
for planets in a volume-limited sample (e.g. Naef et al. 
2010). Twenty-six radial velocities have been collected with 
HARPS with a mean uncertainty of 2.3 ms _I . The dispersion 
of the fit residuals is abnormally large for HARPS (cr _ c = 
3.57 ms -1 ). Because of its faintness, we were not able to 
measure the activity level of HIP 103019. We therefore can- 
not exclude an additional radial-velocity jitter induced by the 
presence of spots on the surface of the star. The fit is satisfac- 
tory (xr = 1-04) if we include additional noise of 3.5 ms" 1 . 
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Fig. 9. Call H emission region for HD 211847 showing re- 
emission at the bottom of the absorption line. 



Table 5. List of targets selected from the literature 



Nr. 


Object 


Nr. 


Object 




HD 


HIP 




HD 


HIP 


1 


13189 


10085 


10 


GJ595 


76901 


2 


30339 


22429 


11 


140913 


77152 


3 


38529 


27253 


12 


162020 


87330 


4 


65430 


39064 


13 


168443 


89844 


5 


91669 


51789 


14 


174457 


92418 


6 


107383 


60202 


15 


180777 


94083 


7 


119445 


66892 


16 


190228 


98714 


8 


131664 


73408 


17 


191760 


99661 


9 


137510 


75535 


18 


202206 


104903 



2.4. Brown-dwarf candidates from the literature 

In addition to the stars from the CORALIE and HARPS surveys, 
we selected 1 8 stars from the literature using the following crite- 
ria: They have announced radial-velocity companions with min- 
imum mass of 13-80 My, were published after 2002, and are not 
flagged as binary stars in the new Hipparcos reduction. The se- 
lected stars are listed in Table 5 and details are given below. 

- HD 1 3 1 89 is a K2 II star with very uncertain mass of 2 -7 M 
and a companion announced by Hatzes et al. (2005) with 
Mosini = 8 -20 My. 

- HD 30339 (F8, 1.1 M Q ), HD 65430 (K0V, 0.78 M Q ), 
HD 140913 (GOV, 0.98 M ), and HD 174457 (F8, 1.19 M ) 
have companions announced by Nidever et al. (2002) with 
Mi sin i ranging between 42 - 78 My. 

- HD 38529 is a G4 IV star with mass 1.48 M Q . Companion 
announced by Fischer et al. (2003) with M 2 sin / = 12.7 My. 
Refined orbits are given by Wright et al. (2009) and Benedict 
et al. (2010). We use the latest orbital elements for analysis. 

- HD 91669 is a K0 dwarf with mass 0.88 ± 0.04 M Q . 
Companion announced by Wittenmyer et al. (2009) with 
M 2 sin / = 30.6 My. 

- HD 107383 is a G8 III giant star with mass 2.7 + 0.3 M Q and 
a companion announced by Liu et al. (2008) with M2 sin i = 
19 AM j. 

- HD 119445 is a G6 III giant star with mass 3.9 ± 0.4 M Q 
and a companion announced by Omiya et al. (2009) with 
M 2 sin / = 37.6 My. 
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- HD 131664 is a G3V star with mass 1.10 + 0.03 M@. 
Companion announced by Moutou et al. (2009) with 
M 2 sin 2= 18.15 My. 

- HD 137510 is a GO IV star with mass 1.42 M and a com- 
panion announced by Endl et al. (2004). The updated data by 
Butler et al. (2006) lists M 2 sin = 22.7 My. 

- GJ 595 is a M3 dwarf with mass 0.28 M . Companion an- 
nounced by Nidever et al. (2002) with M 2 sin i - 60 My. 

- HD 162020 is a K3 V star with mass 0.75 M and a compan- 
ion announced by Udry et al. (2002) with M 2 sin i = 14.4 My. 

- HD 168443 has two companions with minimum masses 
in the planetary (HD 168443b) and brown dwarf domain 
(HD 168443c, M 2 sin / = 18 My) discovered by Marcy et al. 
(2001) and Udry et al. (2004) (see also Udry et al. 2002). 
Updated orbital elements are given by Wright et al. (2009), 
which we use for the astrometric analysis for the outer com- 
panion HD 168443c. 

- HD 180777 is a A9V star with mass 1.7 + 0.1 M and a com- 
panion announced by Galland et al. (2006) with M 2 sin i = 

25 My. 

- HD 190228 is a G5 IV star with mass 0.83 M and a 
planetary-mass companion with M 2 sin i = 3.58 My (Perrier 
et al. 2003). Zucker & Mazeh (2001) performed an astromet- 
ric analysis using Hipparcos and derived M 2 = 67 + 29 My 
at 95 % confidence based on the spectroscopic elements of 
Sivan et al. (2004) (see also Sect. 3.8.2). 

- HD 191760 is a G3 IV star with mass 1.26 + 0.06 M Q . 
Companion announced by Jenkins et al. (2009) with 
M 2 sin i= 38.17 My. 

- HD 202206 is a G6 V star with mass 1.15 M G and has 
two companions with minimum masses in the planetary 
(HD 202206c) and brown dwarf domain (HD 202206b, 
M 2 sin/ = 17.4 My) discovered by Udry et al. (2002) and 
Correia et al. (2005) 

We have now identified our target sample of stars with po- 
tential brown-dwarf companions and have obtained the orbital 
parameters accessible through radial-velocity measurements. As 
the next step, we will search the intermediate astrometric data of 
the new Hipparcos reduction for signatures of the orbital motion 
corresponding to the radial-velocity orbit. 

3. Combined astrometric analysis 

We describe a method to search the Intermediate Astrometric 
Data (IAD) of the new Hipparcos reduction for the orbital signa- 
tures of stars, whose spectroscopic elements are known from a 
reliable radial-velocity solution. The significance of the derived 
orbit is determined through the distribution-free permutation test 
(Good 1994; Zucker & Mazeh 2001). 

Techniques to perform the simultaneous fitting to radial- 
velocity and Hipparcos astrometry data for substellar compan- 
ions and to verify the statistical confidence of the solution have 
been developed by Zucker & Mazeh (2000); Pourbaix (2001); 
Halbwachs et al. (2000); Zucker & Mazeh (2001). More re- 
cently, Reffert & Quirrenbach (2006) and Sozzetti & Desidera 
(2010) used similar approaches to claim the detection of three 
brown-dwarf companions, however with significance not ex- 
ceeding 2-cr. Independently, successful detections of the astro- 
metric signature of exoplanet candidates using the HST fine guid- 
ance sensor (HST FGS) optical interferometer are reported for a 
few stars with spectroscopic elements known from radial veloc- 
ity (e.g. Bean et al. 2007; Martioli et al. 2010). Even though the 
orbital phase is not always fully covered by HST observations, 



for instance the orbital coverage for HD 33636 is 20 % (Bean 
et al. 2007), their superior precision makes the detection possi- 
ble. 

Our analysis is applied to stars with an orbital solution ob- 
tained from radial velocities that provides the spectroscopic el- 
ements P, e, r , a>, and K\. The system is modeled as a pri- 
mary mass Mj orbited by one invisible companion of mass M 2 . 
Using Hipparcos astrometric measurements, we seek to deter- 
mine the two remaining unknown parameters that characterise 
the orbit, which are the inclination i and the longitude of the as- 
cending node Q. The analysis is performed in three steps. First, 
the Hipparcos abscissa, possibly containing the companion sig- 
nature, is constructed from the catalogue data. This is important, 
because the IAD contains only abscissa residuals, which are the 
astrometric residuals after subtraction of the model adopted by 
the new Hipparcos reduction. Second, the best-fit values for i 
and Q. are found by determining the global x 1 minimum on a 
two-dimensional search grid. Finally, the significance of the ob- 
tained solution is evaluated. 

3. 1 . Constructing the Hipparcos abscissa 

Based on the Hipparcos identifier of a given star, the astrometric 
data is retrieved from the new reduction catalogue. It includes 
the astrometric parameters of the solution, the solution type and 
the goodness of fit. The IAD is read from the resrec folder on 
the catalogue DVD of F. van Leeuwen (2007) and contains the 
satellite orbit number, the epoch f, the parallax factor n, the scan 
angle orientation ifr, the abscissa residual <5A, and the abscissa 
error cr A for every satellite scan. 

For the stars in our sample, we encounter the solution types 
'1', '5', '7', and '9'. The solution type '5' indicates that the 
five standard astrometric parameters fit the data reasonably well, 
whereas a solution type ' 1 ' is termed a stochastic solution and 
is adopted when the five-parameter fit is not satisfactory and 
neither orbital nor acceleration models improve the solution in 
terms of x 2 - The solution types '7' and '9' are given, when the 
solution has to include acceleration parameters (proper-motion 
derivatives of first and second order) to obtain a reasonable fit. 
The solution types '1', '7', and '9' can be the first indications 
of a possible astrometric perturbation by an unseen companion. 
The Hipparcos abscissa Amp is then given by: 

AhiP = (a* + Ha* t + A 7 Qf* + Aga*) COS l/f 

+ (6 + us t + A 1 S + AgS) sint/r + mil + SA, ^ 

with the modified right ascension a* = a cos 6, the declina- 
tion 5, the associated proper motions p a * and ps, and the par- 
allax vj. The epoch t is given as time in years since tnip - 
1991.25 A.D. = 2448348.8125 JD and the correction factors for 
solution types '7' and '9' are given by 

A 7 /J=/i p (^-0.8l)/2 (2) 

A 9 p = p p (t 2 - 1.69) r/6, (3) 

where p is either a* or 6 (F. van Leeuwen 2007). fi p and jx p are 
acceleration and change in acceleration in the coordinate direc- 
tion p. For solution types '1' and '5', A-jp and Agp vanish and 
the five parameter model holds. 

Because of the linear nature of Eq. 1, it is usually conve- 
nient to decompose each of the five standard parameters into a 
constant term and an offset, e.g. a* — + Aa* , and to set 
the constant terms to zero to compute the abscissa. As explained 
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in Sect. 3.4, however, our formalism requires the value of the 
parallax m to be present in Eq. 1, whereas we can safely set 
Oq = So = fi a *,o = MftO = for the calculation of the abscissa. 

3.2. Model function 

The model function describes the astrometric signature along the 
scan axis if/ of an isolated star, orbited by one invisible compan- 
ion. To develop our approach, we list the standard formulae (e.g. 
Hilditch 2001) applicable to such a system. The astrometric mo- 
tion of the star is characterised by P, e, To, i, to, £2, and the as- 
trometric semimajor axis a, expressed in angular units. We write 
the Thiele-Innes constants as: 



3.3. Definition of the search grid 

One technique to find a possible orbital signature consists in 
defining a two-dimensional search grid in inclination i and as- 
cending node £2 and to solve for the remaining parameters of the 
model function (Zucker & Mazeh 2000; Reffert & Quirrenbach 
2006). In this way, the goodness-of-fit is obtained for each point 
on the grid and the best solution can be identified as the global 
X 2 minimum. The search grid is defined by i — - 180° and 
Q = - 360°, to avoid negative semimajor axes and sample the 
complete range of values of the cosine-function (cf. Eqs. 4 and 
9). We use a grid size of 30 x 30, which provides a reasonable 
compromise between computation time and grid resolution. 



A — acA c A — cos £2 cos to - sin £2 sin to cos / 

B = acs c B — sin £2 cos to + cos £2 sin w cos / 

F-acf cf = -cos £2 sin w - sin £2 cos w cos / 

G = qcq cq — - sin £2 sin to + cos £2 cos to cos / 



(4) 



where ca,b,f,g are functions of i, £2, and to and a is expressed in 
milli-arcseconds (mas). The elliptical rectangular coordinates X 
and Y are functions of eccentric anomaly E and eccentricity: 



2n 

E-es\nE = — (t-To) 
X = cos E — e 
Y = Vl -e 2 sinE 



(5) 
(6) 
(7) 



If the single-companion assumption is true, the Hipparcos ab- 
scissa is described by a linear equation with 9 parameters, which 
are the 5 standard astrometric parameters plus the 4 Thiele-Innes 
constants: 



Ajy/,p = (a* + fi a * t) cos iff + (6+ fist) sin iff + w II 
+ (BX + G Y) cos iff + (AX + F Y) smif/. 



(8) 



The relationship between astrometrically and spectroscopically 
derived parameters is (Pourbaix 2001): 

a sin / = c P K l P Vl - e 2 cr c P = 3.35729138 ■ 10~ 5 , (9) 

where K\, P, and tv are expressed in ms _1 , years, and mas, re- 
spectively. Hence, the astrometric model function can also be 
written as non-linear function of 12 independent parameters: 

Atf/p = A HIP (a*, 5, w,fi a *,fi s , e, P, T , u>, £2, i, Kij (10) 

The observed radial velocity v rac j of such an orbit is given by 

Vrad = ^1 [C0S(6 I + (jl>) + e cos (l)] + y, (11) 

where y is the systemic velocity. The true anomaly 6 and the 
semi-amplitude K\ are constrained by 



tan - = V(l +e)/(l - e) tan -, 



Ki =(27raisinO/[P(l-e) 1/2 ], 



(12) 



(13) 



3.4. Joint confidence intervals 

For each point on the /-£2-grid we find the least-square solution 
by solving a linear equation using matrix inversion, where e, P, 
To, co, K\ , i, and £2 are fixed. The x 1 value hence obtained for 
each grid-point allows us to derive joint confidence intervals on 
the /-Q-grid. The parallax dependence of the Thiele-Innes con- 
stants is removed by writing 



a = aim, a\ 



and defining the new constant 



cpKxP^ll -e 2 — (14) 
sin/ 



T = a x [(c B X + c c Y) cos if/ + (c A X + c F Y) sin if/] . (15) 

Note that a\ is expressed in AU and corresponds to the angular 
semimajor axis a of the astrometric orbit. Now we can rewrite 
Eq. 8 as 

Amp = (a* + fi a * t) cos if/ + (6 + fi s t) sin^ + m (II + T) . (16) 

This relation is linear in the five remaining free parameters (a*, 
6, Tu, yL a *, fig) and can easily be solved analytically. The com- 
panion signature is solely contained in an additive modulation 
T(z", £2, e, P, To, to, K\) of the parallax factor. It is important to re- 
alise that the presence of T forces the parallax, rather than a par- 
allax offset, to be present in the model function, which originates 
in the abscissa reconstruction from Eq. 1. For the remaining pa- 
rameters it is sufficient to consider the offsets Aa*, A6, Afi a *, 
and Afig to the catalogue values, because the abscissa has been 
constructed with ff ( * = So = fi a *,o = fis,o = 0. 

For every combination of / and £2 the linear equation 16 is 
solved, yielding the five astrometric parameters and the corre- 
sponding^ 2 . The best-fit parameters are identified by the mini- 
mum^- 2 value on the /-£2-grid. The results of the linear adjust- 
ment are not directly used in the quoted final solution, but they 
serve as a consistency check for the results from the non-linear 
fitting and are used for graphical illustration of the joint confi- 
dence intervals as shown in Fig. 12. Low-significance orbits can 
show several local ^-minima corresponding to approximately 
opposite orbit orientations (see also e.g. Zucker & Mazeh 2001; 
Reffert & Quirrenbach 2006). We find that moderate- and high- 
significance orbits have one global ^-minimum and the confi- 
dence contours cover a small area of the /-Q-space (similarly to 
HD 53680 in Fig. 12). The solution parameters corresponding to 
an opposite-orientation orbit are always beyond the 4-cr contour, 
i.e. all orbital parameters of significant solutions are unambigu- 
ously determined. 
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3.5. Determining the final solution 

While the solution method discussed in the previous section is 
favourable because the model function is linear, its accuracy is 
limited by the resolution of the /-£2-grid. For instance, the res- 
olution of a square grid with 900 points is 6° and 12° in i and 
Q, respectively. To avoid this limitation, we determine the best 
solution via ^-minimisation of a model function, where i and 
Q are free parameters and thus can adopt continuous values. 
Therefor, the /-Q-grid is used to define the starting values of a 
non-linear least-squares fit by the Levenberg-Marquardt method. 
The model function is given by Eq. 10, however limited to seven 
free parameters: the 5 astrometric parameters or*, 5, m, jj a *, 
/i,5 plus i and Q. We select the non-linear solution yielding the 
smallest^ 2 and perform 1000 Monte-Carlo simulations. Every 
Monte-Carlo realisation consists of generating a set of Hipparcos 
abscissa measurements and consequent x 1 -minimisation of the 
non-linear model and therefore provides 1000 sets of solution 
elements. 

To incorporate the uncertainties of the spectroscopic param- 
eters we run the complete analysis for 100 sets of spectroscopic 
parameters, where each set is randomly drawn from gaussian 
distributions with mean and standard deviation given by the 
radial-velocity solution and its error, respectively (similarly to 
Sozzetti & Desidera 2010). The final solution accounts for all 
1000 Monte-Carlo solutions, which are obtained for each of the 
100 sets of spectroscopic parameters. To achieve this, we com- 
bine all solutions to yield distributions of 100 000 values for each 
parameter in Eq. 10. The final parameter value and its error is de- 
rived from the mean and the 1-cr confidence interval of the asso- 
ciated distribution, respectively. The values and errors of derived 
quantities, such as the companion mass, are obtained in the same 
way and thus take into account possible correlations between in- 
dividual parameters. 

3.6. Outliers in the astrometric data 

Outlying astrometric measurements are examined if they deviate 
by more than 4-cr from the orbital solution. Here, cr is calculated 
as the root-mean-square of the fit residuals. The datapoint is dis- 
carded if there are at least three measurements at same satellite 
orbit number and if it is more than 3-cr away from these dat- 
apoints. If less than three datapoints at same orbit number are 
present, they are all removed. The analysis is then iterated. This 
procedure is similar to the one applied by Pourbaix & Jorissen 
(2000) and in some cases also agrees with the outlier rejection, 
that is applied by the new Hipparcos reduction itself. Outliers are 
removed for the following objects: GJ 595, HD 3277, HD 43848, 
HD 74842, HD 154697, HD 164427 A, HD 167665, HD 191760, 
and HIP 103019. This never requires more than one iteration. 

3.7. Statistical significance of the orbit 

Orbital solutions can be found for any set of astrometric data and 
the crucial step in the analysis is to determine the credibility of 
the derived orbits. We accomplish this by applying two statistical 
tests: the F-test for its simplicity and the permutation test for its 
lack of underlying assumptions. 

The F-test is extensively used for the statistical analysis 
of astrometric orbit signatures in Hipparcos data (Pourbaix & 
Jorissen 2000; Pourbaix & Arenou 2001; Pourbaix 2001; Reffert 
& Quirrenbach 2006), although it relies on the assumption that 
the measurement errors are Gaussian. That this assumption is not 
necessarily fulfilled for Hipparcos data is mentioned by Zucker 



& Mazeh (2001), who do not apply the F-test for their analy- 
sis. The F-test yields the probability that the null assumption no 
orbital motion present is true, by comparing the ^-values of 
two models with differing number of parameters. In our case, we 
consider the %l -value of the 7-parameter solution by non-linear 
minimisation for each of the 100 draws from the spectroscopic 
elements. The F-test is evaluated on the median ^-value with 
respect to the x\ value, obtained from the standard 5-parameter 
solution, and results in the null probability, which is tabulated 
in Table 8 and gives the probability that the simpler model is 
valid. We use the F-test as an additional indicator for the statis- 
tical confidence of our solution, but the conclusive argument is 
derived from the permutation test. 

The permutation test has the important advantage to be part 
of the distribution-free tests (Zucker & Mazeh 2001) and there- 
fore does not rely on the assumption of Gaussian error distribu- 
tion. For our purpose, it can be summarised by the idea that a 
periodic signal present in a given data set will be destroyed by 
random permutation of the individual datapoints in all but a few 
cases. In the present analysis, the periodicity is defined by the 
radial velocity solution. To derive a significance, the semimajor 
axis of the solution orbit, derived with the original astrometric 
data, is compared to the semimajor axes of pseudo-orbits, which 
are derived from randomly permuted astrometric data. If a sig- 
nal is present in the astrometric data, it will be destroyed by the 
random permutations and the solution orbit will stick out of the 
pseudo-orbits with its large amplitude. 

Because of the pre-reduction present in the original 
Hipparcos IAD (ESA 1997), its permutation had to be done with 
great care (Zucker & Mazeh 2001). Fortunately, it became con- 
siderably easier with the availability of the new Hipparcos re- 
duction IAD. The permutation is applied only to the list of ab- 
scissa residuals <5A, while taking the possible contribution of 
the 7- and 9-parameter solution into account. The remaining 
data, i.e. epochs, parallax factors, and scan angle orientations 
are left unchanged. To obtain the pseudo-orbits, 1000 permu- 
tations are performed and the analysis is run on the synthe- 
sised data, using the nominal spectroscopic elements. The sta- 
tistical significance of the solution orbit is equal to the percent- 
age of pseudo-orbits, that have a smaller semimajor axis than 
the solution orbit (Zucker & Mazeh 2001). Examples for high- 
and low-significance orbits are shown in Fig. 13. In contrast 
to Halbwachs et al. (2000), the analysis of the pseudo-orbits 
is done with spectroscopic parameters identical to the nominal 
ones, hence the method of comparing the distribution of a/o~ a to 
the expected Rayleigh-Rice law is not applicable. 

3.8. Validation of the method 

The combined analysis of radial-velocity and astrometric data 
can lead to erroneous conclusions, if the results are not carefully 
checked (Han et al. 2001; Pourbaix & Arenou 2001). Especially 
the relationship between the inclination and the semimajor axis 
has to be investigated during the interpretation, since over- and 
underestimation of these parameters during the fitting process af- 
fect the derived companion masses. To develop good confidence 
in our results, we validate our method both by simulation and by 
comparison with published work. 

3.8.1. Simulation 

Simulations are employed to investigate the effect of Hipparcos 
measurement precision on the derived orbital parameters and the 
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orbit significance. Of particular interest is the behaviour of the 
occurring biases as function of signal-to-noise given by the ratio 
of orbital semimajor axis and measurement precision. 

Instead of modeling a complete dataset, we chose the case 
of HD 17289, whose orbit is discussed in Sect. 4 and where 
the Hipparcos observations cover twice the orbital period. To in- 
ject the astrometric signature of a companion, we use the satel- 
lite configuration for this star, i.e. the scan angles, the parallax 
factors, and the measurement precision. Noise is generated by 
adding a random bias term to each measurement, which is drawn 
from a Gaussian distribution with standard deviation given by 
the satellite measurement precision. Configurations with differ- 
ent signal-to-noise are generated by altering the radial velocity 
semi-amplitude K\ of the simulated orbit, while keeping all other 
parameters fixed. Consequently, the emulated Hipparcos errors 
in this simulation are Gaussian, which is a simplifying assump- 
tion and may not be fulfilled in reality. However, the simulation 
is valid to investigate the biases of the method. 

The results are illustrated in Figs. 10 and 11, where the ac- 
curacy of the determined values for inclination and semimajor 
axis are displayed as function of statistical significance, derived 
from the permutation test. The accuracy is expressed in terms of 
the difference between determined and true value, i.e. for the in- 
clination it is computed as /derived _ 'true- Generally, the accuracy 
of the derived a and i increase with signal-to-noise and signifi- 
cance. For significance above 2-cr, the accuracy is comparable to 
the precision. However, the semimajor axis retains a small bias 
even at high signal-to-noise, whose magnitude and sign is influ- 
enced by the satellite configuration concerning orbit orientation 
and period. The determination of the orbit inclination greatly im- 
proves with significance. To investigate if our analysis is biased 
even when exceeding the 3-cr significance, we repeated the sim- 
ulation for 3 different satellite configuration, i.e. three different 
stars. We find no sign for a systematic bias, favouring large or 
small inclinations and semimajor axes. Above 3-cr significance, 
derived and true values always agree within the error bars. This 
simulation shows that a large bias can be introduced in the so- 
lution, if an orbit significance of 1-cr is used to qualify a de- 
tection. In contrast, the biases are small if higher significance is 
requested. 

We find that an orbit is detected at better than 3-cr signifi- 
cance, when its semimajor axis amounts to ~70 % of the mea- 
surement precision. This result may depend on the assumption 
of Gaussian errors and the specific satellite configuration, but it 
confirms the warning of Pourbaix (2001), that derived orbits with 
sizes comparable to the instrument precision have to be evalu- 
ated very carefully. 

3.8.2. Comparison sample 

For independent validation of our method and to investigate the 
impact of the new reduction, we screened the literature for ob- 
jects for which a similar analysis was performed. We chose the 
three stars HR 6046, GJ 1069, and HD 190228 with companions 
found by Torres (2007), Halbwachs et al. (2000), and Zucker & 
Mazeh (2001) 3 , respectively, and exhibiting diverse orbit incli- 
nations and companion masses. All these works use the original 
IAD, whereas we use the new reduction IAD. 

For all three stars, good agreement is found between our 
results and the published works, both in terms of value and 
1-cr uncertainty of semimajor axis and companion mass, as is 
shown in Table 6. Also the orbit significances derived from 

3 HD 190228 has also been analysed by Pourbaix & Arenou (2001). 




Significance (o) 

Fig. 10. Simulated difference between the derived semimajor 
axis ^derived and its true value a^. Top: ^derived as function of 
fltme- Filled circles, open circles, and diamonds designate or- 
bits of 1 - 2 cr, 2-3 cr, and > 3 cr significance, respectively. The 
solid line has unity slope and the dashed line shows the single- 
measurement precision of cr A = 7.2 mas. Bottom: The relative 
deviation calculated as (flderived - fltrue)Mrue as function of sig- 
nificance. The accuracy and relative precision of the solution 
improve with significance. At the detection limit of 3-cr, the de- 
viation from the true value is below 0.5 mas. 



pseudo-orbits given by Zucker & Mazeh (2001) for two objects 
match our results. For HD 190228, we find the same solution as 
Zucker & Mazeh (2001) and a 2-cr significance (see Sect. 4.2 for 
further discussion of this object). We confirm the detection of 
the GJ 1069 orbit, caused by a low mass star (Zucker & Mazeh 
2001; Halbwachs et al. 2000), and the binary orbit of HR 6046 
(Torres 2007) at better than 3-cr significance. Especially the case 
of GJ 1069, that is now analysed by three independent teams 
finding essentially identical results, attests that our analysis is 
robust. 

Both the simulation results and the analysis of a sample of 
comparison stars confirm that our method for the combined anal- 
ysis of radial- velocity orbits and Hipparcos astrometry is correct. 
The consideration of the orbit significance ensures that the de- 
rived result is bias-free to reasonable extent. 

3.9. Companion mass limits 

The combination of astrometric and spectroscopic data makes 
it possible to find a complete orbital solution. However, the re- 
sult can be misleading or wrong, if the significance of the solu- 
tion is not verified. For significant orbits the derived parameters, 
e.g. the companion mass, and their errors are valid. For low- 
significance orbits, the derived parameters and their errors are 
very probably false. Therefore, the distribution function for the 
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Table 6. Results for the comparison sample. 



Object 


HIP 




a 

(mas) 






M 2 

(My) 


Significance 

(%) 


Ref. 


HD 190228 


98714 


1 8+ - 5 
1.8_ 05 


[1.8 ± 


0.8] 


65 7+20-4A 

UJ - -20.4 


[67.1--] 


96.3 [95] 


1 


GJ 1069 


19832 


17 7+ o.8 

1 ' ■ ' -0.8 


[17.1 


±0.8] 




[253.5!!^] 


> 99.9 [99.9] 


1 


GJ 1069 


19832 


17 7 +0 - 8 

1 ' ■ ' -0.8 


[17.2 


±0.8] 




[256.7 ± 14.7] 


>99.9 [•■■] 


2 


HD 145849 


79358 


ii r ' 7 


[9.69 


± 0.85] 


1443.5^ 


[1424.8 +7 «] 


>99.9 [•••] 


3 



References. (1) Zucker & Mazeh (2001); (2) Halbwachs et al. (2000); (3) Torres (2007). 

Notes. The values in square brackets are the corresponding values from the reference paper. Our error on the companion mass does not include the 
uncertainty of the primary mass. (o) Primary mass from Sivan et al. (2004) 
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Fig. 11. Simulated difference between the derived inclination ; 
and its true value. The accuracy and the precision of the solu- 
tion improve with significance. At the detection limit of 3-cr, the 
deviation from the true value is below 0.02°. 



companion mass of a low-significance orbit cannot be used to 
derive an upper and lower mass limits. 

The question is then turned to the definition of a significant 
orbit. In the literature, various significance estimators (e.g. F-test 
or permutation test) and detection limits are employed. We use 
a 3-cr detection limit based on the permutation test and deem 
orbits at less than 2-cr unreliable. For a 2-cr orbit we adopt a 
reduced validity and expect that the mass limits do hold. 

In principle, it is possible to derive the upper mass limit for 
any companion by employing simulations similar to the one de- 
scribed in Sect. 3.8.1. Such simulation is extremely elaborate 
and its application to the list of host stars with undetected astro- 
metric motion is beyond the scope of this paper. However, we 
showed that an orbit is detected at better than 3-cr significance, 
when its semimajor axis amounts to ~70 % of the measurement 
precision, if the Hipparcos observations cover one complete or- 
bit. Including an additional confidence margin, we thus claim 
that the signature of a fully covered orbit would have been de- 
tected, if its semimajor axis equalled the single-measurement 
precision. Hence, an upper limit to the companion mass can be 
set by enforcing the equation a — cr A for orbits with significance 
below 2-cr and A^ or b > 1 . 



4. Results 

For all objects in our sample we are able to obtain an orbital 
solution from the combination of the radial-velocity orbit and 
the Hipparcos astrometric data, but not all of these orbits are 
credible. To distinguish a real orbit from meaningless output of 
the fitting procedure, we use the permutation test to evaluate the 
orbit significance. This results in the detection of 6 orbits having 
high significance above 3-cr (Tables 7 and 8). The orbits of 5 
objects have moderate significance between 2-cr and 3-cr. Three 
orbits have low significance between 1-cr and 2-cr and the orbits 
of 17 objects have very-low significance (< 1-cr, Table 10) 4 . 

The Tables 3, 9, and 10 list the spectroscopic and astrometric 
elements of the adopted orbital solution together with the orbit 
significance and give additional information: Af or b is the num- 
ber of orbits covered by the Hipparcos observations, cr A is the 
median Hipparcos single-measurement precision, a sin i is the 
minimum astrometric semimajor axis, M2 (3-cr) are upper and 
lower companion-mass limits at 3-cr, a K \ is the semimajor axis 
of the relative orbit, re( ] is me reduced chi-square value of the 
adopted 7-parameter solution, and the Null probability gives the 
result of the F-test. The quoted errors correspond to Monte-Carlo 
based 1 -cr-confidence intervals. 

For orbits with significance below 2-cr and provided that the 
orbital period is fully covered by Hipparcos observations, we 
derived the upper limit to the companions mass M2, up -iim by en- 
forcing the equation a — cr A (cf. Sect. 3.9). 

4. 1 . Orbits with high significance exceeding 3-cr 

Clear orbital signatures are detected for six target stars. The or- 
bit significances exceed 99.7%, which we adopt as criterion for 
undoubtful detection. The orbits show a large variety in size 
(a = 2.3 - 35 mas) and inclination (i = 12 - 173°) and re- 
veal companion masses ranging from 90 My to 0.52 M . The 
obtained relative precision on the companion mass is typically 
of the order of 10 % and does not include any contribution of 
the primary-mass uncertainty. The visualisation of Hipparcos as- 
trometric data is not very intuitive, because positions are mea- 
sured in one dimension only, i.e. along the scan-angle orienta- 
tion 1//. We choose a representation similar to Torres (2007) and 
the resulting stellar astrometric orbits are shown in Fig. 20. For 
better display, we compute normal points for each satellite or- 
bit number. Their values and errors are given by the mean and 
standard-deviation of the Hipparcos abscissae obtained during 
one given satellite orbit, respectively. Normal points are used 

4 As motivated in Sect. 2.3, we do not present an astrometric analysis 
for HD 4747 and HD 21 1847 and therefore the number of objects listed 
here adds up to 31 instead of 33. 
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Fig. 12. Joint confidence contours on the 
i-£l-ghd for the low-significance orbit of 
HD 167665 (left) and the high-significance or- 
bit of HD 53680 (right). The contour lines 
correspond to confidences at 1-cr (solid), 2-cr 
(dashed), 3-cr (dotted), and 4-cr (dash-dotted) 
level. The crosses indicate the position of the 
best non-linear adjustment solution for each of 
the 100 draws and the star corresponds to the 
adopted orbit. 
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Fig. 13. Histogram showing the semimajor 
axes of the 1000 pseudo-orbits for the low- 
significance orbit of HD 167665 (left) and the 
high-significance orbit of HD 53680 (right). 
The vertical solid line indicates the semimajor 
axis of the non-permuted, best-fit solution and 
the dotted line indicates the median Hipparcos 
single-measurement precision <r A for the re- 
spective star. 



Table 7. Astrometric solution elements for targets with orbits that have significance above 1-cr. 



Object 


Aa* 


AS 




Act 


Afi a * 


A^ 


i 


Q. 




(mas) 


(mas) 


(mas) 


(mas) 


(mas yr" 1 ) 


(mas yr -1 ) 


(deg) 


(deg) 


HD 17289 


-6 7 +a8 

"•'-0.8 


6 +a8 

"-0.8 


20 l +u - b 


-0.1 


, 9+O.6 
1 - Z -0.6 


9 9+0.7 
Z - Z -0.7 


173.2^- 4 


2 5+ 1 " 


HD 30501 


9 r 30 


-12 6 +35 

u - -3.5 


48 l +0 - 4 


0.2 


-3 +0 - 8 

-'■"-0.8 


-4 +18 

^•"-1.8 


49. i:^ 1 


no 7+12.0 
"•'-12.4 


HD 43848° 


2 7+1.3 

^•'-1.3 


4 0+ 1 - 3 

-1.3 


26.4+°; 7 


-0.1 


-0 2 +1 ° 
"■^--l.o 


-5 6 +1 1 


165 -°-lo 


97Q 9+12.9 


HD 43848* 


3 7 +L9 

-'•'-1.9 


7 Q+2.0 
' - y -2.0 


26.5^7 


0.1 


-0.7^3 


-4.4!!;! 


155.7!« 4 6 7 


237.9 + _\\ 1 1 


HD 53680 


6 2+ 2 - 1 

-2.1 


41 3 +4 '° 


58.2™ 


0.8 


-21.5+£* 


2.0!!; 7 


163.6![- 4 


238.9! 2 ; 9 


HD 164427 A 


4 +0 - 7 

u - H -0.7 


1 5 +0 - 5 


25.8** 


-0.5 


o r 06 


2 +0 4 

"•^-0.4 


li.s:!! 


338.9!! 2 ; 3 


HIP 103019 


5 6 +L6 

-1.6 


4 2+1.3 

^•^-1.3 


30.2-[;l 


1.9 




2 +16 

"• -1.6 


160.92? 


143.1^« 


HD 3277 


-1 o +0 - 5 


3 +0 ' 5 

"-'-0.5 


riA q+0.6 


0.1 


"0-8!^ 


-0 6 +0 5 

"•"-0.5 


166.5! 2 -; 


270.2+!! 2 


HD 131664 


-0 3 +L6 

"-1.6 


1 Q+0.8 
i -^-0.8 


18.6^ 


0.9 


-1 6 +1 1 
i -"-i.i 


9 q+0.9 
z - J -0.9 


167.7^ 


321.0^ 


HD 154697 


4 2 +L7 


5 9 +2 - 8 


?q 4+1-0 


-1.1 


7 7+1.8 
'•'-1.9 


9 -5+1.8 

^•-'-l.s 


148.8+ 4 ;* 


293.9^ 


HD 174457 


-0 7 +0J 

"• ' -0.7 


i 9+0.6 
i --'-0.6 


18-5^7 


-0.2 


-0 3 +0 8 

—0.9 


-0-8^ 


137 +11 - 9 

"'■"-14.1 


133.0!! 72 


HD 190228 


6 +0 - 6 

"•"-0.6 


-1 3 +0 - 7 


16 +0 " 


-0.2 






4 3 +L8 

^•■'-1.0 


61.0™ 


HD 74014 


-5 7 4+ 33 - 6 

-"-33.5 


AQ 4+39.6 
-39.3 


29 5 +0 - 7 


-0.0 


-13.333 


_4 Q+5.2 


171 l +3 ° 

1 ' 1 - 1 -10.6 


248. 7^ 7 ; 7 


HD 168443c 


-0 7 +0 " 7 
"•'-0.7 


5 +0 - 6 

"•-'-0.6 


26 3 +0 - 7 


-0.4 


1 2 +L0 

1 - -0.9 


-o r 10 


35.7!i 4 | 


130 1+33.4 


HD 191760 


-0 5 +L0 
"• J -i.o 


1 o +0 - 8 

1 -"-0.8 


11 5 +L1 


-0.8 


7 +L4 

"•'-1.4 


_0 4+1.0 

"•^-1.0 




Q+46.6 



Notes. The top 6 targets have high-significance orbits (> 99.7 %). The bottom 3 targets have low-significance orbits (1-cr- 2-cr) and the parameters 
represent the formal solution. The remaining 5 orbits have moderate significance (2-cr -3-cr). {a) Final solution at lower eccentricity. {V) Formal 
solution at higher eccentricity. 
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Table 8. Solution parameters and companion masses for stars with orbits that have significance above 1-cr. 



Ohiert 


Sol. 


orb 


IT . 

u A 


a sin i 






2 \~ ) 




Y 1 

/(7,red 


Null nrnh 


Si tmificancp 




type 




(mas) 


(mas) 


(mas) 


(Mj) 


(Mj) 


(mas) 




(%) 


(%) 


HD 17289 


1 


2.1 


7.2 


1.29 


10.8IJ" 


547.4^ 


( 386, 700) 


32.0 


0.32 


7e-29 


> 99.9 


HD 30501 


5 


0.6 


3.3 


10.44 


i4.o:;i 


89.6!| 2 ; 3 


(65, 132) 


147.7 


0.79 


0.005 


99.9 


HD 43848" 


9 


0.5 


4.2 


2.32 


8-9!!;; 


ioi.8*j|5 


(59, 150) 


90.8 


1.09 


4e-8 


> 99.9 


HD 43848* 


9 


0.5 


4.2 


4.93 


10.51U 


120.2!}^ 


(73, 179) 


91.7 


1.07 


le-8 


> 99.9 


HD 53680 


1 


0.7 


5.7 


9.73 


34.9«; 2 


226.7!£° 


( 159, 301) 


160.7 


0.82 


3e-16 


> 99.9 


HD 164427A 


5 


11.1 


2.4 


0.49 


9 O+0.5 


269.9!«; 2 


( 101,482) 


12.8 


1.15 


0.003 


99.7 


HIP 103019 


5 


1.2 


7.1 


3.33 


io.s + _\i 


188.1^ 


(118,280) 


53.4 


1.13 


4e-10 


> 99.9 


HD 3277 


5 


24.8 


3.6 


0.58 


2 5 +0 - 5 


344.2™ 


( 143, 600) 


9.4 


1.00 


5e-4 


97.7 


HD 131664 


5 


0.6 


4.7 


0.88 


A -X+d-t 
^-0.8 


89.611^ 


(44, 153) 


60.3 


0.77 


2e-5 


97.9 


HD 154697 


7 


0.6 


3.1 


6.26 


11 7 +L7 


151 9 +249 


( 81,235) 


89.2 


1.60 


le-6 


96.0 


HD 174457 


5 


1.3 


3.2 


1.77 


2-8^6 


107.8! 23 ;* 


(65, 190) 


35.2 


1.22 


12.3 


95.9 


HD 190228 


5 


1.0 


3.5 


0.13 




4Q 4+14.8 


(10,99) 


32.9 


1.08 


0.3 


95.4 


HD 74014 


5 


0.1 


2.5 


9.84 


64 +32 3 

ut.u _ 334 


491 o+243.8 
-252.1 


(51,1568) 


235.4 


1.34 


19.9 


80.6 


HD 168443 


5 


0.6 


1.9 


1.25 


9 4+0.6 


34.3!»;° 


(18,67) 


75.9 


1.53 


3.6 


75.0 


HD 191760 


5 


1.6 


3.3 


0.46 


1 Q+0.8 


185.3*" 


(39, 659) 


16.0 


1.39 


27.6 


77.5 



Notes. The top 6 orbits have high significance > 99.7 % and are therefore fully characterised. The bottom 3 orbits are of low significance (l-er-2-cr) 
and represent the formal solution. The remaining 5 orbits have moderate significance (2-<x -3-cr) and provide lower and upper mass-limits for the 
companions. (fl) Final solution at lower eccentricity. {b) Formal solution at higher eccentricity. 

Table 9. Spectroscopic elements of literature targets with orbits of better than 1-cr significance. 



Object 


HIP 


Mi 


M 2 sin i 


P 




e 


Ki 




T 








Ref. 






(M ) 


(Mj) 


(day) 






(ms~ 


') 


(ID*) 






(deg) 




HD 131664 


73408 


1.1 


18.2 


1951.0 ±41.0 


0.64 


±0.02 


359.5 ± 


22.3 


52060.0 ± 


41.0 


149.7 


± 1.0 


1 


HD 174457 


92418 


1.2 


65.8 


840.8 + 0.1 


0.23 


±0.01 


1250.0 4 


: 10.0 


52020.0 ± 


:4.0 


139.0 


± 1.0 


2 


HD 190228 


98714 


0.8 


3.6 


1146.0+16.0 


0.50 


±0.04 


91.0 ± 


5.0 


51236.0 ± 


25.0 


100.7 


±3.2 


3 


HD 168443c 


89844 


1.0 


18.1 


1748.2+ 1.0 


0.212 + 


: 0.001 


298.1 i 


:0.6 


53769.8 ± 


0.0 


64.7 


±0.5 


4 


HD 191760 


99661 


1.2 


38.0 


505.6 ± 0.4 


0.63 


±0.01 


1047.8 i 


:38.7 


54835.7 ± 


2.1 


200.4 


±0.3 


5 



References. (1) Minniti et al. (2009); (2) Nidever et al. (2002); (3) Perrier et al. (2003); (4) Wright et al. (2009); (5) Jenkins et al. (2009). 



only for the visualisation and not during the orbit adjustment. We 
prove that orbit detections are also possible in apparently diffi- 
cult conditions, such as incomplete orbit coverage by the satellite 
(N or b = 0.5 for HD 43848) and an instrument precision compa- 
rable to the the orbit size (HD 164427 A). The final 7-parameter 
fit for the detected orbits has an acceptable reduced chi-square 
value (Xj Kd = 0.8 - 1.2) in all but one case. For HD 17289 this 
value is very small at 0.3 and the Hipparcos precision is unusu- 
ally large at 7.2 mas, which may indicate that the astrometric 
errors for this star are overestimated. As for the radial-velocity 
measurements, the light emitted by the actual stellar companion 
of this star may have disturbed the Hipparcos observations. 



Three out of these six stars have non-standard Hipparcos so- 
lutions. Thus, first indications for the astrometric disturbance in- 
duced by the companion were already detected during the stan- 
dard Hipparcos data analysis. An additional outcome of the com- 
panion solution are refined positions, distances, and proper mo- 
tions of the objects (Table 7). For 5 stars the parallax estimate 
becomes more precise, whereas its value is compatible within 
the error bars. Only the parallax of HIP 103019 is 1.9 mas larger 
than the value given by the new reduction. Changes in position 
and proper motion are of the order of 1 mas and 1 mas yr~' , with 
extrema of 41 mas and 22 mas yr~', respectively. 



- HD 17289: The orbit inclination is 173.2±0.5° and the com- 
panion mass is 0.52 ± 0.05 M , which makes it the most 
massive companion detected in this study. HD17289 is thus 
identified as a binary star with mass ratio q — 2. Goldin & 
Makarov (2007) used solely Hipparcos astrometry to find an 
astrometric orbit with semimajor axis a = 1 1 .7 + j which is 
compatible with our more precise value of a — 10.8 ± 0.7. 
Because the orbit is completely characterised, we can calcu- 
late the projected angular separation between the two com- 
ponents. Along the orbit, their separation on the sky ranges 
between 15-49 mas. The mass-luminosity relations of 
Delfosse et al. (2000) estimates an apparent V-magnitude of 
the secondary of V% = 13.18. The primary magnitude ob- 
tained from Hipparcos is V = 7.56 (Table 1) and the ex- 
pected intensity contrast between the secondary and primary 
is ~1:180 in the visible. The combination of close separation 
and moderate intensity contrast causes a small and periodic 
bias in the radial-velocity measurement, see Sect. 2.3. 

- HD 30501: The orbit inclination is 49*g° deg and the com- 
panion mass is 90 + 12 My, which is the lightest companion 
detected in this study. This companion is just at the threshold 
between a very low-mass star and a brown dwarf. 

- HD 43848: The orbit inclination is 165^ deg and the com- 
panion is a very low mass star with mass 0.10 ± 0.01 M . 
This value is slightly lower and more precise than the 
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Table 10. Parameters for targets with very-low significance orbits (<l-cr) 



Object 


HIP 


Sol. 


N„rb 


OA 


Mi 


Mi sin i 


a sin ( 


Ref. 


Null prob. 


Significance 


^^2,up-lim 






type 




(mas) 


(M ) 


\Mj) 


(mas) 




(%) 


(%) 


(Me) 


HD 52756 


33736 


5 


22.5 


4.2 


0.8 


59.0 


0.53 


1 


26.1 


54.0 


0.72 


HD 89707 


50671 


5 


3.5 


2.7 


1.0 


53.0 


1.54 


1 


19.5 


52.0 


0.11 


HD 167665 


89620 


5 


0.2 


1.8 


1.1 


50.6 


7.70 


1 


25.1 


54.0 




HD 189310 


99634 


5 


82.0 


4.2 


0.8 


26.0 


0.10 


1 


2.6 


58.3 


2.93 


GJ595 


76901 


1 


12.2 


19.9 


0.3 


60.0 


3.50 


2 


89.5 


8.0 


0.63 


HD 13189 


10085 


5 


1.9 


3.5 


4.5 


14.0 


0.01 


3 


71.3 


11.0 


9.62 


HD 30339 


22429 


5 


60.0 


3.8 


1.1 


77.8 


0.12 


2 


35.6 


23.0 


10.46 


HD 38529 


27253 


5 


0.4 


1.6 


1.5 


13.4 


0.80 


4 


55.9 


9.0 




HD 65430 


39064 


5 


0.3 


3.1 


0.8 


67.8 


12.79 


2 


21.0 


61.0 




HD 91669 


51789 


5 


2.1 


6.0 


0.9 


30.6 


0.44 


5 


17.9 


52.0 


0.53 


HD 107383 


60202 


5 


3.5 


1.1 


2.7 


19.4 


0.10 


6 


64.2 


9.0 


0.22 


HD 119445 


66892 


5 


2.7 


2.2 


3.9 


37.6 


0.06 


7 


77.3 


0.0 


1.86 


HD 137510 


75535 


5 


1.4 


1.9 


1.4 


22.7 


0.70 


8 


45.7 


34.0 


0.061 


HD 140913 


77152 


5 


7.3 


5.1 


1.0 


43.2 


0.49 


2 


20.7 


20.0 


0.56 


HD 162020 


87330 


5 


92.7 


5.0 


0.8 


14.4 


0.05 


9 


20.8 


50.0 


7.32 


HD 180777 


94083 


5 


40.8 


1.2 


1.7 


25.0 


0.11 


10 


51.9 


7.0 


0.27 


HD 202206b 


104903 


5 


3.9 


3.8 


1.1 


17.4 


0.26 


11 


48.7 


10.0 


0.27 



References. (1) This work; (2) Nidever et al. (2002); (3) Hatzes et al. (2005); (4) Benedict et al. (2010); (5) Wittenmyer et al. (2009); (6) Liu et al. 
(2008); (7) Omiya et al. (2009); (8) Butler et al. (2006); (9) Udry et al. (2002); (10) Galland et al. (2006); (11) Correia et al. (2005). 



one derived by Sozzetti & Desidera (2010), who obtain 
O.II^q q^Mq based on the radial velocities of Minniti et al. 
2009. This solution is found with the lower-eccentricity or- 
bit (cf. Sect. 2.2). For completeness, we include the for- 
mal astrometric solution for the less probable and higher- 
eccentricity radial-velocity orbit in Tables 7 and 8. The re- 
sulting companion mass obtained in this case would be 
higher at M2 — 0. 1 1 + 0.02 M . Hipparcos astrometry cannot 
be used to determine the actual eccentricity, because the fit 
quality in terms of \ 2 is comparable in both cases. 

- HD 53680: The orbit inclination is 163.6^ deg and the 
companion is a low mass star with mass 0.22 ± 0.02 M Q . 

- HD 164427A: The orbit inclination is 11.8+|] deg and the 
companion is a low mass star with mass 0.26 ± 0.06 M . 
The stellar companion has also been detected by Zucker & 
Mazeh (2001), who derived a slightly less inclined orbit and 
therefore obtained the larger mass of 0.35 + 0.09 M G . 

- HIP 103019: The orbit inclination is 160.9^ 7 deg and the 
companion is a very low mass star with mass 0. 18 ± 0.03 M Q . 

4.2. Orbits with moderate significance between 2-cr and 3-cr 

The orbits of 5 stars are found with significance of 95.4-99.7 %. 
Although we do not consider that these solutions accurately 
characterise the stellar orbits, the moderate level of significance 
indicates that the orbital motion is present in the astrometric 
data. The formal orbits are shown in Fig. 21 and allow us to ob- 
tain a visual impression of the fit quality. Four out of five orbits 
have sizes comparable to or smaller than the instrument preci- 
sion and the astrometric signal is therefore hidden. The permu- 
tation test allows us to assess the credibility of the derived orbit 
especially in this difficult regime. Since there is signal in the as- 
trometric data, the orbit inclination is constrained and we expect 
the companion mass to be confined within the 3-cr limits given 
in Table 8. 



- HD 3277: The companion is of stellar nature with its mass 
confined within 0.14-0.57 M . The formal companion mass 
at moderate significance is 0.33 + 0.07 M . The upper mass- 



limit which we independently derived from a — cr A is 
A^2,up-iim = 0.5 1 M Q and compatible with our solution. 

- HD 131664: The companion is a very low-mass star or 
a brown dwarf with a mass of 44 - 153 My. The avail- 
able data is not sufficient to distinguish between these two 
cases. The formal companion mass at moderate significance 
is 90± 16 My. Although the orbit significance agrees with the 
result of Sozzetti & Desidera (2010), who obtain a compan- 
ion mass of 23^ 6 My with 95 % confidence, the mass range 
derived from our analysis is significantly higher. 

- HD 154697: The companion is a low-mass star with a mass 
confined within 0.08 - 0.22 M . The formal companion mass 
at moderate significance is 0.15 + 0.02 M Q . The longitude 
of the ascending node Q of this orbit is hardly constrained 
and is correlated with the position offsets Aa* and A6. The 
Q- Ac* -correlation is shown in Fig. 14. Because the adopted 
solution is derived from the mean of the complete distribu- 
tion, it does not represent a valid solution with a valid orbit, 
as is also visible in Fig. 21, where the fit residuals are large. 
However, this does not affect the validity of the derived com- 
panion mass, because that is independent of £1 

- HD 174457: The companion is a very low-mass star or a 
brown dwarf. Its mass is confined within 0.06 - 0.18 M . 
The formal companion mass at moderate significance is 
0.10 ± 0.02 M Q . The upper mass-limit which we indepen- 
dently derived from a = <x A is M2, up -ii m = 0.12M Q and com- 
patible with our solution. 

- HD 190228: The companion is probably a brown dwarf, al- 
though we cannot exclude a very low-mass star. Assuming 
a primary mass of 0.83 M (Perrier et al. 2003), we find the 
formal companion mass of 49 ± 15 My with lower and upper 
limits of 10 - 99 My, cf. Table 6. The upper mass-limit de- 
rived from a = cr A independently reproduces this result and 
yields M 2 , up -i im = 99 My. 



4.3. Orbits with low significance between 1-cr and 2-cr 

The orbits of three stars are found at low significance of 68.3 - 
95.4 %, i.e. the astrometric data contains very little or no or- 
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Fig. 14. Correlation between the adjusted position and longi- 
tude Q present in the distribution of the combined 100000 
Monte-Carlo realisations for HD 154697. The adopted solution 
is marked with the diamond. 



bital signal. As we have shown by simulation (cf. Sect. 3.8.1), 
orbital solutions at this significance level are prone to large bi- 
ases. Therefore, the solution parameters shall not be considered 
valid in these cases. The formal solution of the fitting procedure 
is given for completeness in Tables 7 and 8. 

- HD 74014: The poor Hipparcos coverage of 10 % of this 
long-period orbit impedes the detection of the astrometric 
orbit. With a minimum semimajor axis of a sin i = 9.8 mas, 
the orbit would have been detected on a longer timespan, 
given the instrument precision of 2.5 mas. Because the orbit 
is not fully covered, we cannot set an upper companion-mass 
limit M 2 , U p-iim- 

- HD 168443: The formal solution for the companion mass 
of 34 ± 9 My is in agreement with the result of Reffert & 
Quirrenbach (2006), who derive 34 + 12 My, using the radial 
velocities of Marcy et al. (2001) and the original Hipparcos 
data. From our analysis, using the CORALIE radial velocities 
and the new Hipparcos reduction IAD, we conclude that the 
mass obtained for the outer companion HD 168443c is of 
low confidence. Because the orbit is not fully covered, we 
cannot set an upper companion-mass limit M2. up -iim. 

- HD 191760: The minimum orbit size of this star is very small 
(a sin i = 0.46 mas) and Hipparcos astrometry does not re- 
veal the signature of the companion discovered by Jenkins 
et al. (2009). However, we can set an upper mass-limit for 
the companion of M2, up -iim = 0.28 M . 

4.4. Orbits with very-low significance and the detection of 
one brown-dwarf companion 

The analysis of 17 objects yields orbits with significance below 
1-0". Because the astrometric data does not contain orbital in- 
formation, the derived result is physically meaningless and only 
a mathematical solution. Therefore we do not list the solution 
parameters. Table 10 shows the targets and their basic parame- 
ters. Two characteristics can be observed that impede the astro- 
metric orbit detection. In the one case, applicable to HD 65430 
and HD 167665, the orbital coverage is poor and the orbit is 
not detected although the expected signal is large compared 
to the instrument precision. In the other case, for instance for 



GJ 595, HD 13189, and HD 140913, the minimum orbital sig- 
nature given by a sin i is very small compared to the instrument 
precision. 

The companions of six stars have maximum mass M2 iUp _ii m 
above 0.63 Ms, which does not represent a considerable con- 
straint of the system. Nine companions have upper mass-limits 
of M2 ;Up _ii m = 0.11- 0.63 M and thus in the M-dwarf mass 
range. Finally, the mass of the companion of HD 137510 is con- 
fined between M 2 sin/ = 22.7 My and M2, up -iim = 64.4My. 
Therefore, we find that HD 137510b has to be a brown dwarf. 
This detection became possible, because of the confidence we 
gained on the capability of Hipparcos to detect astrometric or- 
bital motion. Using a simpler argument, Endl et al. (2004) de- 
rived an upper mass-limit of 94 My, which was not stringent 
enough to prove the substellar nature of the companion. 

This sample also contains HD 38529, for which Benedict 
et al. (2010) found the orbit inclination of the outer companion 
HD 38529c and derived its mass of M 2 = 17.6 My using HST 
FGS astrometry. Our analysis made use of the radial velocity or- 
bit given by Benedict et al. (2010), but fails to detect a signifi- 
cant orbit, although the Hipparcos precision of 1.6 mas is in the 
range of the HST orbit size (a = 1.05 ± 0.06 mas). However, the 
Hipparcos orbit coverage is poor at 0.4, which can explain the 
non-detection. 



5. Discussion 

The recent compilation of close companions to Sun-like stars 
with minimum masses of 10 - 80 My by Sozzetti & Desidera 
(2010) lists 39 objects. With the discovery of 9 new compan- 
ions, we increase the known number of such objects by about 
20 %. Ten companions in the Sozzetti & Desidera (2010) list 
are probably stars. We confirm the stellar nature in two cases 
(HD 43848B and HD 164427B) and detect 6 new stellar com- 
panions (HD 3277B, HD 17289B, HD 30501B, HD 53680B, 
HD 154697B, and HIP 103019B) at the 3-cr significance level. 
All of these companions are M dwarfs with masses between 
90 My and 0.52 M . For 3 companions, we derive mass limits 
that enclose the boundary between low-mass stars and brown 
dwarfs: HD 131664 (M 2 = 44 - 153 My), HD 174457 (M 2 = 
65 - 190 My), and HD 190228 (M 2 = 10 - 99 My). HD 190228 
initially was a planet-host star discovered by Sivan et al. (2004). 
We analysed it to validate our method and consequently recog- 
nised its importance as a promising brown dwarf candidate. 

Furthermore, we set upper mass limits below 0.63 M Q to the 
companions of nine stars. These companions are therefore M- 
dwarfs or brown dwarfs. The upper mass-limits of 12 compan- 
ions within our sample remain unconstrained, either because the 
derived limit is very high (> 0.63 M ) or because the orbital 
period is not covered by Hipparcos measurements. Finally, we 
found that the companion of HD 137510 is a brown-dwarf with 
a mass range of 22.7 - 64.4 My. 

Our results are summarised in Fig. 15, where we show the 
companion-mass constraints obtained from Hipparcos astrome- 
try for all targets in our analysis sample. 

5. 1 . Combining radial velocities and Hipparcos astrometry 

We used the astrometric data of the new Hipparcos reduction to 
set constraints on the masses of potential brown-dwarf compan- 
ions, detected in radial-velocity surveys. Radial-velocity mea- 
surements of a stellar orbit yield 5 of the 7 Keplerian elements 
and do not constrain the orbit's inclination and longitude of the 
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Fig. 15. Constraints on the companion masses of the target sam- 
ple. Open squares show the minimum companion mass M2 sin i 
derived from radial velocities. Blue filled circles indicate the for- 
mal companion masses M2 derived from astrometry for orbits 
with >2cr-significance. Grey-shaded areas show the companion 
mass interval constrained by astrometry. The vertical dashed line 
indicates the adopted threshold mass of 80 Mj between brown 
dwarfs and stars. Twelve systems could not be constrained by 
Hipparcos astrometry, either because the derived mass limit was 
very high (>0.63 M ) or because the orbital period was not cov- 
ered by Hipparcos measurements. 



ascending node. We demonstrated how these two remaining ele- 
ments can be determined, if the orbital signature is present in the 
Hipparcos astrometry, and how to distinguish a significant orbit. 

We used the permutation test to define the criterion for orbit 
detection. The implementation of this powerful test is consider- 
ably simplified by the use of the new Hipparcos reduction, be- 
cause the individual measurements of the IAD are independent. 
A source of possible errors and complications that was present in 
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Fig. 16. Comparison of the orbit significance indicators derived 
from the permutation test and from the F-test. 



the original Hipparcos data is therefore avoided. We emphasize 
the importance to eliminate outliers of the IAD, because other- 
wise the final result can be falsified. On a couple of comparison 
targets drawn from the literature, we have verified that our anal- 
ysis method is robust. On these targets, the performance of the 
new and the original Hipparcos IAD is comparable and no major 
difference of the results was observed. 

In addition to the permutation test, we simultaneously per- 
formed the F-test and can thus compare the performance of both 
significance indicators, keeping in mind that the F-test is valid 
under the assumption of Gaussian errors only. Figure 16 dis- 
plays the respective values listed in Tables 8 and 10. We find 
a general agreement between the two indicators, i.e. the null 
probability decreases with increasing permutation significance. 
However, the permutation test is more stringent. If we had cho- 
sen to rely on the F-test and imposed a null probability of 0.3 
% as detection limit, we would have detected 10 orbits. These 
overlap with the 6 orbits, which we detected with the criterion 
based on the permutation test (cf. Table 8). 

A problem of similar studies in the past was the appearance 
of biases towards small orbit inclinations and thus towards large 
companion masses (Han et al. 2001; Pourbaix 2001). Using sim- 
ulations, we showed that our method is not affected by such bias 
when using an appropriate detection limit in terms of orbit sig- 
nificance. Figure 17 shows the distributions of orbital inclination 
and semimajor axis for all orbits with better than 2-cr signifi- 
cance. No preference of very small inclinations can be observed, 
as it would be expected from biased solutions (Pourbaix 2001). 
However, we note that no orbit is detected at inclination above 
50°. This is explicable, because in this study we have selected 
stars with minimum companion masses of M2 sin i = 13 -80 My. 
As shown by Grether & Lineweaver (2006), the frequency of 
brown-dwarf companions is very low compared to the stellar 
companions and therefore small inclinations are more likely to 
be detected. 

We showed that our analysis is ultimately limited by the as- 
trometric precision of the satellite and therefore cannot constrain 
orbits with a size below typically a = 2 mas, which in most 
cases will impede the detection of a true brown dwarf compan- 
ion. In addition, the Hipparcos mission duration of about 1000 
days limits the periods of detectable orbits to ~2000 days. 

In summary, we have developed a new method to combine a 
radial velocity solution with Hipparcos astrometric data that has 
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Fig. 17. Inclinations and astrometric semimajor axes of the or- 
bital solutions presented in Tables 7 and 8. The inclination is 
displayed as angular deviation from a face-on orbit (correspond- 
ing to i = 0° or i = 180°). Filled symbols indicate the orbits with 
high-significance above 3-cr, while open circles stand for orbits 
with moderate significance between 2-cr and 3-cr. 



several advantages compared to previous similar approaches: it 
uses the new Hipparcos reduction, which simplifies the math- 
ematical formulation and the use of the permutation test com- 
pared to the implementation of e.g. Zucker & Mazeh (2001) with 
the original Hipparcos data. The new method defines rigorous 
detection limits in terms of orbit significance and yields com- 
panion mass limits for low-significance orbits, provided they are 
covered by astrometric measurements. It is stringent enough to 
avoid the output of biased results (an example is the case of 
HD 38529, see Reffert & Quirrenbach 2006 and Benedict et al. 
2010). Because of the general format of the new reduction IAD, 
our method is universal and applicable to any astrometric data 5 . 

5.2. Brown-dwarf companions in the CORALIE survey 

The CORALIE planet-search sample contains about 1600 Sun- 
like stars within 50 pc. After a duration of 12 years, 21 potential 
brown companions have been observed in this survey. Eighteen 
are discussed in this work and three additional hosts of poten- 
tial brown-dwarf companions are listed by Sozzetti & Desidera 
(2010) 6 and contained in the CORALIE sample. Assuming that 
most candidates contained in this survey have already been 
found, we can derive a typical frequency of 1.3 % for poten- 
tial brown-dwarf candidates in close orbit (< 10 AU) around 
Sun-like stars. Seven of these candidates are identified as stellar 
companions through this study. Additionally, Zucker & Mazeh 
(2001) showed that the companions of HD 18445, HD 112758, 
and HD 217580 are of stellar nature. Thus, at least 10 out of 21 
brown-dwarf candidates are in fact low-mass stars, which cor- 
responds to a rate of 48 %. Based on the uniform stellar sample 
surveyed by the CORALIE planet search, we thus obtain an upper 



5 This is true provided that the implicit approximations are valid, i.e. 
the astrometric precision is larger than about 1 mas, see e.g. Konacki 
et al. (2002). 

6 Three more stars from this list were not included because their min- 
imum mass is M2 sin i < 13 Mj. 
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Fig. 18. Cumulative mass distribution of potential brown-dwarf 
companions in the CORALIE survey. The blue dashed line shows 
the distribution of all 21 candidates. The black solid line shows 
the distribution for the 1 1 remaining candidates after removal of 
the 10 stellar companions. The companion of HD 38529 with 
M2 - 17.6 My (Benedict et al. 2010) is included. For compar- 
ison, the red dash-dotted line shows the expected cumulative 
distribution if all companions had masses of 80 My under the 
assumption that the orbits are randomly oriented in space. 



limit of 0.6 % for the frequency of close-in brown-dwarf com- 
panions around Sun-like stars within 50 pc, which is close to 
the frequency of < 0.5 % that was initially quoted by Marcy & 
Butler (2000). 

Figure 18 shows the distribution of A^sim in the brown- 
dwarf mass range for the CORALIE sample. While the cumula- 
tive distribution (blue dashed line) of all 21 candidates is rea- 
sonably compatible with a linear increase, i.e. a flat distribu- 
tion function, the situation changes significantly after removal 
of the 10 stellar companions. As illustrated in Fig. 19, most stel- 
lar companions have Mi sin/ > 45 My, and after their removal, 
a cumulative distribution with particular shape emerges (black 
solid line in Fig, 18). It exhibits a steep ascent in the range of 
13 - 25 My, which contains one half of the objects, followed by 
a slower increase up to 60 My. Surprisingly, there is no compan- 
ion left with minimum mass larger than 60 My. 

In summary, the emerging cumulative distribution does not 
support a uniform distribution function of companion masses. 
The possibility, that all companions in the 13 - 25 My range are 
in fact stellar companions seen at low inclination, is both sta- 
tistically unlikely and also unplausible because their astrometric 
signature would be large and probably detectable. Instead, these 
companions supposedly have M2 ~ M2 sin i and may represent 
the high-mass tail of the planetary distribution. 

We note that the distribution found here, which is based on 
radial velocities, is compatible with the first brown-dwarf com- 
panions found in transiting systems, such as CoRoT-3b (M2 = 
21.7 My, Deleuil et al. 2008) and CoRoT-15b (M 2 ~ 60 My, 
Bouchy et al., in preparation). 
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Fig. 19. Histogram of minimum masses of potential brown- 
dwarf companions in the CORALIE survey. The light-grey his- 
togram shows the distribution of all 21 candidates. The dark-grey 
histogram shows the distribution of the 1 1 remaining candidates 
after removal of the 10 stellar companions. 



5.3. Metallicities of the brown-dwarf host stars 

We showed that the companion of HD 137510 is a brown dwarf. 
Additionally, the companion of HD 190228 can be considered a 
brown dwarf at 95 % confidence. The host stars have metallic- 
ities of [Fe/H] = 0.373 dex (Butler et al. 2006) and [Fe/H] = 
-0.24 dex (Perrier et al. 2003), respectively. Compared to the 
mean metallicity <[Fe/H]>= -0.09 dex with root-mean-square 
°"<[Fe/H]> = 0.24 dex of the uniform sample of 451 nearby Sun- 
like stars from Sousa et al. (2008), the two brown-dwarf hosts are 
therefore metal-rich in one case and unremarkable in the other 
case. The question whether hosts of brown dwarf companions 
show a particular metallicity distribution, as the planet hosts do 
(Santos et al. 2001), cannot be answered with this small sample 
and thus remains open. 

6. Conclusions 

We presented radial velocity solutions of 15 stars with potential 
brown dwarf companions, of which 9 are new discoveries. At 
this point, the CORALIE and HARPS surveys have increased the 
number of known potential brown-dwarf companions of Sun- 
like stars by 20 %. 

Eight companions in our search sample were identified as 
being of stellar nature with masses above 80 My. These include 
one companion with a mass of 90 My and thus very close to the 
brown-dwarf mass boundary. The upper and lower mass limit for 
three companions was derived, but include the star-brown dwarf 
mass boundary and can therefore not be used to identify the ob- 
ject's nature. Yet, the companion of HD 190228 is a brown dwarf 
at 95 % confidence. Upper mass limits were set for the compan- 
ions of 9 stars and revealed a brown dwarf orbiting HD 137510. 
The masses of 12 companions could not be constrained. 

Our findings are in accordance with the presence of the 
brown-dwarf desert, because only 3 brown dwarfs are identi- 
fied in our search sample of 33 companions: HD 190228b (at 
95 % confidence), HD 137510b, and HD 38529b (Benedict et al. 



2010). In contrast, 8 brown-dwarf candidates were identified as 
M dwarfs. Based on the CORALIE planet-search sample, we ob- 
tain an upper limit of 0.6 % for the frequency of brown-dwarf 
companions around Sun-like stars and confirm their pronounced 
paucity. We find that the companion-mass distribution function 
is rising at the lower end of the brown-dwarf mass range, sug- 
gesting that in fact we are seeing the high-mass tail of the plan- 
etary distribution function. 

We find good agreement with the previous similar studies 
of Halbwachs et al. (2000) and Zucker & Mazeh (2001). Our 
results are affected by the Hipparcos astrometric precision and 
mission duration, which limits the minimum detectable mass of 
the combined analysis. There are indications, that some of the 
remaining potential brown-dwarf companions are actual brown 
dwarfs. Those could be identified in the future, when astrome- 
try at higher precision will be available. An example is the case 
of HD 38529c, whose mass we were not able to constrain with 
Hipparcos, but that was identified as brown dwarf by Benedict 
et al. (2010) using HST astrometry. 

Higher-precision astrometry will be required to obtain an 
accurate census of low-mass companions in close orbit around 
Sun-like stars and to derive the relative frequencies of planets, 
brown dwarfs, and low-mass stars. The next major step in this 
direction is the GAIA astrometry satellite, which is the successor 
of Hipparcos and will make it possible to determine the masses 
of many radial-velocity companions. 
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Fig. 20. Visualisation of the high-significance orbits. Top panels: Astrometric stellar orbits projected on the sky. North is up and 
East is left. The solid red line shows the orbital solution and open circles mark the individual Hipparcos measurements. Dashed 
lines with orientation along the scan angle </> and length given by the O-C residual of the orbital solution connect the measurements 
with the predicted location from our model. The blue solid circles show the normal points for each satellite orbit number. The curl 
at the lower left corner indicates the orientation of orbital motion. Bottom panels: O-C residuals for the normal points of the orbital 
solution (filled blue circles) and of the standard 5-parameter model without companion (black crosses). The error bars of the normal 
points correspond to the dispersion of Hipparcos measurements if there are several per satellite orbit and to the individual Hipparcos 
abscissa error if there is only one measurement. 
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Fig. 21. Visualisation of the formal solutions for the moderate-significance orbits. The panels are analogous to Fig. 20. The formal 
orbit of HD 154697 is not valid and exhibits large residuals (see the text for explanation). It is shown here to illustrate the possible 
problems that can occur when visualising astrometric orbits. Despite their moderate significance, the four remaining orbits visually 
appear well constrained. 
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